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1. Introduction 

Nuclear physics is a well-established and mature field of study. A large number of experiments 
have been performed during the past decades that have led to a great understanding of nuclei and 
nuclear processes. We are currently moving into an exciting time for the field, one in which, for 
the first time, properties of nuclei will be rigorously determined directly from QCD. This will first 
occur through numerical calculations of quantities using lattice QCD, matching onto effective field 
theories (EFTs) that have been developed during the last decade or so, and using the EFT's to 
compute the properties of the light nuclei. The first lattice QCD calculation of the deuteron will 
be a milestone for nuclear physics. It is this scenario that has motivated a number of us to leave 
the relative comfort of our "analytic" offices, in which we developed the EFT's and fretted over the 
ever increasing number of counterterms, and to attempt calculations of multi-nucleon systems with 
lattice QCD. 

2. Motivations 

There are a number of compelling reasons to use lattice QCD to perform nuclear physics 
calculations. While it is important to explore and understand the building-blocks of nuclei, the 
proton and neutron, a greater challenge lies in understanding the structure and interactions of nuclei. 

2.1 The Emergence of Nuclear Physics from QCD 

One of the deep issues in our field, one for which we currently have little or no understanding, 
is how nuclear physics emerges from QCD. The Lagrangian describing QCD is simple: there is a 
gauge-covariant kinetic term for the quarks, a mass-term for the quarks and the Yang-Mills term 
describing the glue. There are a small number of fundamental constants that appear in this La- 
grangian: Aqcd and the six quark masses. This Lagrangian and these seven constants (along with 
the electroweak interactions) give rise to all of nuclear physics. We wish to establish a rigorous 
pathway from QCD to nuclei. As mentioned previously, this will be accomplished by perform- 
ing lattice QCD calculations of some quantities, determining counterterms in the appropriate low- 
energy EFT from these calculations, and then using the EFT to compute quantities of interest in 
relatively simple systems. This EFT will be subsequently matched onto the many-body machinery 
that has been developed and used by nuclear physicists, which will then be used to compute the 
properties of more complex systems 1 . 

As a simple example, consider the classic nuclear process np — > dy. This is a process in 
which a low-energy neutron is captured by a proton to form a deuteron, emitting M 1 radiation in 
the process. As the energy increases from threshold the El amplitude rapidly increases and soon 
dominates the Ml amplitude. A cartoon of the process at the quark-diagram level is shown in 
fig. [T]. The Maiani-Testa theorem ^ prevents direct "computation of this process" on the lattice. 
Instead, one must compute other quantities on the lattice and use these quantities to obtain the 
matrix elements of interest. This will be accomplished by matching a lattice calculation to an 
EFT calculation. For the process at hand, one can either use the pionless theory, in which only 

'This may involve a lattice calculation of the effective field theory itself, e.g. Ref. ||l|^[^[0|, without reference to 
QCD. 
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Figure 1: A cartoon of the QCD calculation of the simplest nuclear radiative capture process np — > dy. 

nucleons and photons are dynamical degrees of freedom, or the pionful theory, in which nucleons, 
pions and photons are dynamical, see Ref. [0J. The pionless theory is very useful for processes at 
momentum much below the pion mass, and this theory has been explored out to several orders in 
the power-counting. The pionful theory is somewhat less developed, but if one is interested in the 
quark-mass dependence, as will be the case for the initial calculations performed at quark-masses 
greater than those of nature, the pionful theory will be used. In either theory, there are contributions 
that are unrelated to the NN scattering amplitude, described by gauge-invariant operators involving 
nucleons and the electromagnetic field strength tensor. The coefficients of the operator(s) will be 
determined by matching to a lattice calculation. For this process {np — > dy) there are experimental 
constraints that provide a very precise determination of these counterterms, and a lattice calculation 
will provide (other than the quark mass dependence) nothing more than a check of the method. 

Ultimately, moving toward more complex systems one will need to match these EFTs, the zero, 
one and higher nucleon EFTs to the many-body techniques that have recently been developed. The 
nuclear shell-model continues to move toward a rigorous underpinning based on the renormaliza- 
tion group. In principle, one starts with a very large model space, as defined by the number of 
levels of the harmonic oscillator (HO) (for instance) and matches chiral perturbation theory and 
NN EFT's to observables in this model space. One then systematically reduces the model space, 
by integrating out the high-levels of the HO. This is exactly analogous to the renormalization group 
methods we are used to dealing with in particle physics, except this is now constructed for bound 
states. As states are integrated out, the Hamiltonian and operators are renormalized. Eventually 
one is left with a small, low-energy space -the shell-model space-and a large excluded space. The 
operators and Hamiltonian evolve in such a way that the energy-levels and matrix elements in the 
shell-model space reproduce those computed in the full-space, see Ref. [^]. 

2.2 Calculations of Processes where Experiments are not Possible 

Lattice QCD will have significant impact on nuclear physics by computing quantities for which 
there is little or no data. One such example of this is to perform calculations that impact our 
understanding of dense nuclear matter. I do not mean asymptotically dense, but at the densities that 
arise in the collapse of stars, that are a few times nuclear matter density. A second example is the 
calculation of matrix elements of operators that probe physics beyond the standard model, such as 
j3/3 -decay. 
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2.2.1 Supernova, Neutron Stars and Black Holes 

The nature of the remnant of supernova 1987A is still undetermined. It is unknown if the 
remnant is a neutron star or a black hole (and it is unknown if it underwent spherical collapse or 
aspherical collapse). At the heart of the matter is nuclear physics, and in particular the compress- 
ibility of nuclear matter at densities that are a few times that of ordinary nuclear matter. To resolve 
this issue, one needs to know the composition of nuclear matter in this regime, and this is where 
lattice QCD can contribute. 

As one compresses nuclear matter composed of neutrons, protons and electrons, the Fermi- 
energy of the electronic component increases rapidly compared to that of the nucleons. At some 
point it becomes energetically favorable for the system to increase the density of neutrons via the 
weak interaction, rather than to further increase the density of electrons and protons, despite the 
proton-neutron mass-splitting. Hence, the name neutron star. Thus it is the neutrons and their 
interactions that primarily dictate the compressibility of nuclear matter at densities higher than 
nuclear matter densities. This discussion ignores the possibility that other hadronic components 
can become energetically favored. In 1986 Kaplan and Nelson []|, ^] pointed out that the K~n 
interaction should significantly reduce the mass of the K in neutron matter, and as stressed by 
Bethe and Brown [10], if the K mass falls below the electron Fermi energy, the formation of a 
condensate of K '~'s, (K~) ^ will occur. Nuclear matter would then be composed of neutrons, 
protons and a kaon condensate which has a softer equation of state than pure neutron matter. The 
nK~ scattering amplitude is poorly known, and thus there is considerable uncertainty in the density 
dependence of the K mass. This ambiguity can be resolved by a careful lattice QCD study. 

A competing process arises from the density dependence of the £~ mass. If the mass of the 
£~ drops below the electron chemical potential, then nuclear matter will be composed of neutrons, 
protons and £~'s. For obvious reasons, this system will have a softer equation of state than pure 
neutron matter. The nL~ scattering amplitude is also poorly known, only accessible through the 
properties of hypernuclei, and so one can presently only speculate as to the existence of such matter. 

Lattice QCD can have significant impact on determining the equation of state of hadronic mat- 
ter at high densities. The two processes discussed above, kaon-condensation and hyperonic matter 
can be refined by a calculation of the neutron-kaon and neutron-hyperon scattering amplitudes. 



2.2.2 j3j3 -Decay 

The last decade has seen remarkable progress in our understanding of the neutrino sector. The 
measurements by SuperK and SNO demonstrating that neutrinos have non-zero mass, and that 
the mass eigenstates are not the flavor eigenstates, confirmed what had been speculated for many 
years based upon the results of the Davies chlorine experiment. It remains to be determined if the 
neutrinos are Dirac or Majorana. In the later case, lepton number is not conserved, as the neutrino 
is its own anti-particle. Observation of a lepton number violating process, such as neutrinoless 
/3/3 -decay, would be yet another sensational discovery in the neutrino sector. 

Nature has provided a limited number of nuclei that undergo /3/3 -decay, and the constraints 
on building detectors to observe j3/3 -decay allow for only a handful of elements to be considered. 
Perhaps the most widely known is Germanium, where one looks for a peak in the e~e~ invariant 
mass spectrum resulting from 76 Ge — > 76 Se e~e~, sitting on top of the background from the lepton 
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number conserving process 76 Ge ~^ 16 Se e~V e e-V e as a signal of lepton-number violation. As 
the simplest baryonic process that can contribute to this process is of the form pp — ► nne~e~, and 
higher body operators are expected to be suppressed, it is clear that a rigorous calculation of the 
rate of such decay is not simple. There are complications at the many-body level, where the calcu- 
lation of the nuclear rate from the two-nucleon operator (local or non-local) is notoriously difficult. 
Further, for a massive neutrino, or where the process is dominated at short-distances, matching 
from QCD onto the nuclear EFT suffers from issues similar to those that arise in the calculation 
of nonleptonic weak matrix elements [11 , 12]. Lattice QCD is the only way to rigorously explore 
these matrix elements. 



2.3 The Dependence of Nuclei and Nuclear Processes on the fundamental constants of 
Nature 

Before we can be content with our understanding of nuclei and nuclear processes we must 
know and understand how they depend upon the fundamental constants of nature. The recent 
suggestion that the electromagnetic coupling is time dependent, based upon the splitting between 
atomic absorption-lines in gas clouds at different red-shifts (along the line of sight to quasars) [pj|], 
renewed interest in such a possibility. I look forward to an independent verification of this exciting 
result. 

If we consider a nuclear process that is of central importance to the development of life, e.g. 
3a ^ n C, we should be able to fully understand how the rate for this process, and the rate for the 
production of carbon in stars depends upon the mass of the light-quarks, the scale of the strong 
interaction and the value of the electromagnetic coupling constant. At this point in time we have 
very little control over any aspect of this calculation. One of the reasons for this is that this process 
is highly fine-tuned, but also we have only a limited knowledge about how nuclear physics depends 
upon the quark masses (even in the absence of fme-tunings). Clearly it is not going to be easy to 
map out the spectrum of 12 C as a function of the quark masses, however, this is a worthy goal for 
nuclear physics. 



2.4 Fine-Tunings in Nuclear Physics 

During the last few years the role of fine-tunings in physics has become a central focus of 
particle physics. The question of why the cosmological constant, A ~ 10~ 47 GeV , is so small and 
yet nonzero has sparked great interest. However, fine-tunings also play a central role in the strong 
interactions. 

The nucleon-nucleon potential can be dissected into roughly three length scales. The long- 
distance part is well described by one pion exchange (OPE). Performing a phase shift analysis of 
NN scattering data while treating the mass of the pion as a free parameter, along with its couplings, 
recovers the physical pion mass and its axial coupling. What is more impressive is that a recent 



analysis of NN scattering [ ]14[ ] has shown that the two-pion exchange contribution computed in chi- 
ral perturbation theory provides a better fit than the traditional "a-meson" exchange. What's really 
impressive is that the values of the local counterterms that enter at this order are consistent with 
those extracted from the single-nucleon sector. Therefore both the long-distance and intermediate- 
distance component of the NN-potential is well described by chiral perturbation theory. While 
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many "old-school" NN potentials have a one meson exchange short-distance repulsion, the modern 
potentials have a short-distance component that is optimized in shape and strength, consistent with 
EFT power-counting, to best fit the data. 

The truly fascinating aspect of this, is that while the NN potential has three quite distinct 
distance scales, the s-wave wavefunction resulting from solving the Schrodinger equation for the 
two-nucleon system is approximately horizontal at the edge of the potential, i.e. the scattering 
lengths are unnaturally large, ~ —24 fm in the 'Sb-channel and ~ +5.4 fm in the Si -channel (the 
channel containing the deuteron, which is bound by ~ 2.2 MeV). There is clearly a fine-tuning 
among the regions of the potential in order for this to be true. Clearly, this is all controlled by 
the parameters of QCD, but in nature these parameters have conspired to place us very near an 
infrared fixed-point in the hadronic sector (if the scattering lengths were infinite the system would 
exhibit scale-invariance and would be at an infrared, unstable fixed-point) 1 . Proximity to this fixed 
point means that the EFT used to describe such systems does not follow the usual power-counting 
developed for systems that have natural scattering lengths, where it is simply a matter of counting 
engineering dimensions. The four-nucleon operators contribute at the same order in the counting 
as the chiral limit of OPE, while deviations from the chiral limit are suppressed [jT6|]. In order to 
have states near threshold, the coefficients of the four-nucleon operators are fine-tuned to achieve 
the delicate cancellation between kinetic and potential energy that exists. 

The 3a -> 12 C process that was discussed previously, is also finely-tuned. In order for this pro- 
cess to proceed with any appreciable rate there are two important fine-tunings that have occurred. 
First, the ground state of % Be is barely unbound, and in the stellar interior one has 2a Be, and 
the s Be state lives for an unnaturally long period of time, enough time for a third a to interact to 
form 12 C through radiative capture. However, the rate would still be insufficient except for the fact 
that there is a level in 12 C at precisely the right energy to significantly enhance the capture rate. It 
is quite easy for me to believe that we can only understand this set of fine-tunings via the anthropic 
principle, as without carbon we would not be having this discussion. However, I have been unable 
to reconcile these fine-tunings with those in the NN-sector, as I have no reason to believe that the 
form of the NN interaction required to greatly enhance the triple-alpha process would also put us 
close to a fixed-point in the NN-sector. There has been no work to try to reconcile these features, 
and it is something that should be understood. 



3. The Quark- Mass Dependence of Nucleon-Nucleon Scattering 

The development of EFTs for nuclear physics during the last 15 years since Weinberg's pio- 
neering papers on the subject, have allowed us to rigorously explore the quark-mass dependence 
of simple nuclear systems. We have learned that simple analytic dependences on the quark masses 
are incorrect and the actual quark mass dependence is quite complex. This is compounded by the 
fine-tunings discussed in the previous section. 

During the last couple of years, the quark-mass dependence of the two-nucleon sector has been 
explored ][l7|, |l8|, 19]. The scattering lengths as a function of the pion mass are shown in fig. ||. 



At the first non-trivial order, in addition to the quark mass dependence from OPE, there are also 



2 The fact that both channels, and not just one, are close to the fixed point is understood in the \axge-N c limit of 
QCD [|l^] leaving us to understand why the system as a whole is near the fixed-point 
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Figure 2: The NN scattering lengths in the s-wave channels as a function of the pion mass at NLO in the 
NN EFT. The two regions defined by the points (green and black) correspond to two different choices of the 
(sets of) coefficients of the m^-dependent four-nucleon operators jl8|], The three (blue) solid squares are the 
scattering lengths resulting from a quenched lattice QCD calculation [|20||. 



contributions from four-nucleon operators that have a single insertion of the quark mass matrix with 
coefficients, the D2's, that are currently unknown. To determine the reasonable ranges of scattering 
length at a given quark mass dimensional analysis has been used to constrain the coefficients of 
these operators. Fig |2| suggests that it is likely that the di-neutron remains unbound as one moves 
toward the chiral limit, while it could be bound or unbound as the mass is increased. In the deuteron 
channel, the deuteron binding energy moves toward its natural value for decreasing quark mass, but 
again, it could be bound or unbound at higher masses. In fig. ^| I have shown the only lattice QCD 
calculation of this system [20]. This pioneering lattice calculation is quenched and at large pion 
masses, outside the range of the EFT's that currently exist, and therefore, unfortunately, does not 
shed light on the question at hand. 

The program to undertake is to perform a lattice QCD calculation of the scattering lengths (or 
of the deuteron) at the larger masses that are currently accessible, but inside the range of validity of 
the EFT, match to the and use the EFT to make the chiral extrapolation to map out the entire 
region of interest. 



4. Two Hadrons on the Lattice 



Motivated by the physics described in the previous sections two collaborations have emerged: 
NPLQCD and StellaQCD, see fig. || 3 . The current focus of these two efforts is to extract the scat- 
tering lengths, and more generally, the phase shifts associated with low-energy nucleon-nucleon, 
nucleon-hyperon and hyperon-hyperon scattering. 

4.1 The Maiani-Testa Theorem (The End of the Innocence) and Luscher's Method 

The first hurdle that one encounters when considering the extraction of a scattering ampli- 

3 The participating institutions are University of Washington, University of New Hampshire, Lawrence Berkeley 
Laboratories, College of William and Mary, University of Groningen and University of Barcelona. 
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Figure 3: The logos of NPLQCD and StellaQCD. 



tude from a lattice QCD calculation is the Maiani-Testa theorem. Essentially no-one outside of 
the conference attendees appreciates the implications of this theorem, and typically the impression 
that most have is that one will be able to calculate everything directly from lattice QCD. Unfortu- 
nately, this is false. The Maiani-Testa theorem states that one cannot extract two-body (and higher) 
S-matrix elements from Euclidean space Green functions at infinite-volume except at kinematic 
thresholds. Obviously, this has huge implications for the computation of nuclear processes, even 
the simple ones like np — > dy, with lattice QCD. 

In the two-body sector, Liischer was able to arrive at a method to circumvent this constraint [ 21 , 
22]. Below inelastic threshold the volume dependence of the two-particle energy levels is deter- 
mined by the scattering amplitude. Thus determining the energy-levels of two particles in a finite- 
volume allows for a determination of the scattering amplitude (at the location of the energy-levels). 
Liischer's formula relating the scattering amplitude to the location of two-particle states at finite- 
volume is 



pcotd(p) 



where 



1 

%L 




(4.1) 



s(r?; 



4k A 



j ' 



(4.2) 



and it is understood that the UV cut-off, Ay — > oo. In particular limiting cases one can arrive at 
analytic expressions for both the continuum states and bound states. For a/L — ► 0, the lowest 
energy state is located at (we are using the nuclear physics definition of the scattering length) 



Eq = + 



Ana 
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(4.3) 



where the coefficients are c\ = —2.837297, C2 = +6.375183, and the bound state, if one exists, is 
located at 
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where (pcotd)' = pcotd evaluated at p 2 = -y 1 . The quantity 7 is the solution of 



7 + pcotS|_2 = _ 







(4.5) 



In the case of a/L — ► °° (more precisely Lpcotd — > 0), of relevance to the NN system, one finds 
that 



An 1 
ML? 



[d\ + d% Lpcotdo + ... ] , 



(4.6) 



where the coefficients are d\ = —0.095901, = +0.0253716 and where /?cotc>o is evaluated at an 
energy E = ^d\. 

4.2 The Present and the Future for NN Scattering 

One of the exciting aspects of studying two-nucleons on the lattice is that rigorous calculations 
can be performed today, at unphysical pion masses, even if the systems have scattering lengths that 
are much larger than the lattices. The reason for this is that it is the range of the interaction that 
is relevant, and not the scattering length. The range of the NN potential is dictated by the mass of 
the pion, and while the physical pion mass requires volumes larger than > 5 fm, a pion of mass 
~ 350 MeV requires volumes larger than > 2.5 fm, which means that meaningful calculations can 
be performed on the available MILC lattices. 



2 

q m 



small-L expansion 
— ■ Luscher approximation 
% exact solution 



50 



L (fm) 



— bound-stale approximation 
— - Luscher approximation 

— — small-L expansion 
C exact solution 




Figure 4: The position of energy levels of two nucleons in a cubically-symmettic finite volume of spatial 
dimension L at the physical value of the pion mass. The y-axis is q 2 ~ E„L 2 . The points correspond to exact 
solution of eq. (4.1), while some of the curves correspond to the approximations in eq. (4.3), eq. (|4.4|), and 
eq. (4.6). As the interactions between nucleons in the two channels are similar, the spectta in small volumes 



are also similar. Further, as the scattering lengths are large in both channels (a fine-tuning between potential 
and kinetic energy), the spectra in large volumes are similar, except for the bound state in the channel with 
positive scattering length. 



The simplest tool for analyzing lattice results is the pionless EFT, which is appropriate below 
the pion inelastic threshold. Beane, Bedaque, Parreno and I [^] explored what would be found 
in a lattice QCD calculation of the NN system at finite volume at the physical value of the pion 
mass using the well known NN scattering amplitude. The results are shown in fig. |j. The energy 
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of levels associated with two nucleons in the [ Sq and the mixed 3 5i — 3 Di channels were computed, 



and solutions to eq. ( |4.1| ) (containing both the continuum states and bound state) are shown by 
the solid circles in fig. ||. As expected the bound state solution asymptotes at large volume to the 
deuteron binding energy, while the lowest continuing levels move toward vanishing energy 

To explore NN interactions in smaller volumes, one needs to solve the pionful theory in a 
finite-volume with periodic boundary conditions (or whatever boundary conditions one chooses, 
for instance, twisted boundary conditions as first suggested by Bedaque fl24| ] and explored further 
in Refs. [ p5 , 26]). This has not yet been performed, and given that the pionful NN EFT calculations 
are performed to lower orders than the pionless calculations, this is going to require some serious 
efforts by the NN EFT community. 



4.3 Electroweak Interactions 

Perhaps the simplest non-trivial process in nuclear physics involving electroweak interactions 
is the radiative capture process np — > dy. This is a classic nuclear physics process, and provides 
a clean demonstration of meson-exchange-currents (MEC's) (electroweak interactions that are not 
related by gauge symmetry to the NN scattering amplitude). Lattice QCD is going to be able to 
compute this process by matching onto the low-energy EFT (either pionless or pionful), determin- 
ing the coefficients of the gauge-invariant operators that are not related to the NN scattering ampli- 
tude, and using the EFT to compute the amplitudes of interest. Detmold and I explored what could 
be learned from lattice QCD calculations of processes such as this using background-field tech- 



niques [27]. We found an expression for the energy-levels of two nucleons in a finite volume in the 
presence of a background magnetic field, and also a background Z° field (for the weak-dissociation 
of the deuteron). The background fields mix the y Sq and the 3 S\ states, and the spectrum is some- 
what complicated. However, there is at least one level in the spectrum (and the challenge is to find 
lattices for which it is the ground state or a low lying excitation) that is sensitive to the coefficient 
of the gauge-invariant operator. This sensitivity comes about because the deuteron and the near 
bound state are finely-tuned, and placing the system in a finite volume disturbs this fine-tuning. 
The contribution of the local operator can restore the fine-tuning in the presence of the background 
field. It is in this scenario that the level becomes sensitive to the value of the coefficient. 

For np — > dy, the energy of two nucleons in a background field (in a cubically symmetric 
lattice are) 



where 
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(4.8) 



The Ui 



contributions result from the one-nucleon interactions with the background field, while the two- 
nucleon-background field interactions are described by l\. It turns out to be desirable to work 
with asymmetric volumes to suppress the Landau-level structure, and fig. ^ shows the levels for 
particular values of the parameters. For np — > dZ *, the relation becomes 
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Figure 5: NN energy levels in an asymmetric volume in the presence of a background magnetic field (left 
panel) and background Z°-field (right panel) [p7[]. There is at least one level that is sensitive to the value of 
the gauge-invariant operator involving two nucleons and the background field. 
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5j = S(p 2 + wj) , S 2 = S{p 2 - 
The energy-levels in an asymmetric volume are shown in fig. |5| 
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4.4 The AqAq Potential 

An observable that may impact our understanding of nuclear forces is the potential between 
two Aq baryons in the heavy-quark limit. The potential between to B-mesons in the heavy-quark 
limit has been explored previously [28, 29, [3l| ], both in quenched and dynamical calculations. 
This interaction also resembles the nuclear force, including OPE. For two Aq baryons, isospin 
symmetry forbids contributions from OPE, and the leading long-distance behavior arises from the 
exchange of two pions (TPE). While TPE leads to a better description of the NN scattering phase- 
shift data than heavy meson exchange, one would like to better understand this component of the 
interaction. The TPE contribution to the Aq potential is only part of the TPE contribution to the 
NN potential, as the Weinberg-Tomazawa term is absent, the axial matrix elements are different, 
and the spectrum of the Aq — Eg sector differs from the N — A sector. However, one will be able 
to directly explore what TPE "looks like" in this system, for instance, if this part of TPE well- 



described by perturbation theory. Arndt, Beane and I computed this potential J32[ ] in #PT and also 
in partially-quenched #PT. As expected it "falls like a rock", as the relevant mass-scale is 2m K , a 
plot of which can be found in fig. |6[ 



5. Nonleptonic Hyperon Decays 

As the nonleptonic decays A — > /?7T , Z~ — > nn~, £ + — > nn + , £~ — ► An~ and their isospin 
partners have been well studied for an extensive period of time, the experimental data and theoret- 
ical efforts to understand the data can be found in many textbooks. Recently we emphasized that 
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Figure 6: Logarithm of the potential between two Ag's. 
this is a problem worthy of exploration with lattice QCD and examined the amplitudes using SU(2) 



symmetry [33] 



While the axial matrix elements g& ~ 1.27 and g\z ~ 0.60 are well-known experimentally, 
only an upper limit currently exists for Z~ — ► iPeV, and thus there is no direct measurement of the 
coupling gzz, which contributes to the £-decay p-wave amplitudes. Orginos has been awarded time 
by SciDAC to compute this matrix element, and the other octet-baryon axial matrix elements. 

Inserting numerical values for the couplings and masses into the S— wave and P— wave ampli- 
tudes gives the numerical values for the nonleptonic amplitudes shown in table ||. At leading order, 

Table 1: Weak Amplitudes in 51/(2) jPT at LO (gzz = 0.30 -> 0.55) 



Decay 


Theory 


g/W Expt 


srf^ Theory 


i/( p ) Expt 


A — > p7t~ 


1.42 (input) 


1.42 ±0.01 


0.56 


0.52 ±0.02 


£~ — * n%~ 


1.88 (input) 


1.88±0.01 


-0.50 -0.14 


-0.06 ±0.01 


L + — > H7t + 


0.0 


0.06 ±0.01 


±0.42 -» ±0.08 


1.81 ±0.01 



the P— wave amplitude in A decays is well predicted, as it is in the three-flavor theory. However, 
we do not expect significant modifications to this result from higher orders in the two-flavor theory. 
Further, we see that the P— wave amplitude for £~ — > n%~ is quite sensitive to the value of g^x., 
which is presently quite uncertain. At the upper limit of the allowed range for gzz, the P— wave 
amplitude is close to what is observed. Finally, we see that the P— wave amplitude for £ + — ► n7t + 
is not close to the experimental value for any reasonable value of gzz, and theory underestimates 
the experimental value by ~ 4 in the best case. 



6. NPLQCD 

The NPLQCD collaboration applied for, and was granted, time from SciDAC for exploratory 
investigations of the AW, NY and YY systems (y=hyperon), and also to extract the a-term and 
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strong isospin breaking in the nucleon mass. We were awarded ~ 8% of JLab resources, ~ 
40 Gfiop — yrs to perform computations with domain-wall-ferrnions[|4|, 35] on the available MILC 



lattices, using Chroma/QDP++ developed by Edwards and his team at JLab [|36|]. We have ana- 
lyzed I = 2 Tin scattering and have computed the correlation functions for AA. In what follows we 
present the results obtained by performing contractions with DWF propagators generated by the 
LHPC collaboration, who very kindly allowed us to use them in our work. 

6.1 1 = 2 TTTT-Scattering 

Before tackling the two baryon systems we thought it wise to explore 1 = 2 7T7r-scattering [j37|]. 
Not only is this a good "warm-up" exercise, but it is also interesting physics in itself. 
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Figure 7: Logarithm of the ratio of the 1 = 2 %n correlation function to the square of the % correlation 
function projected onto the momentum |p| = state (left panel) and the |p| = 1 state (right panel). For 
the excited state (right panel) we have shown the function with and without subtracting the ground state 
contribution. 

The correlation functions shown in fig. were used to extract the ground state and first ex- 
cited state energies of on 20 3 x 64 staggered MILC lattices, with a physical dimension 
~ 2.5 fm and lattice spacing ~ 0.125 fm. The physical pion masses from these calculations are 
m n ~ 290, 350, and 480 MeV. 

In fig. H we show the dimensionless quantity m^a 2 versus the dimensionless quantity m n /f % . 
The lowest mass CP-PACS point computed with Wilson fermions is also shown. It is remarkable 
how well the tree-level result from chiral perturbation theory describes the lattice calculation. Also, 
shown is the residual after removing the tree-level amplitude and multiplying by ~ f^/rn^. The 
lattice data is not precise enough to observe the anticipated chiral logs that enter at next order, 
however, at these quark masses the counterterm and the chiral log essentially cancel. Clearly, a 
more precise study of this system is required in order to identify the chiral logs that are predicted 
to be present, and to allow for a significantly more precise chiral extrapolation to compare with 
experiment. Our extrapolated value is 



m n a2 



-0.0426 ±0.0006 ±0.0003 ±0.0018 , 



(6.1) 
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Figure 8: The results of our lattice QCD calculation of m n a% as a function of m n /f K (ovals) |]37[]. Also shown 
are the experimental value from Ref. p8[ ] (diamond) and the lowest quark mass result of the dynamical 
calculation of CP-PACS [[39| (square). The gray band corresponds to a weighted fit to our three data points 
using the one-loop J-PT formula. 



where the first error is statistical, the second is the systematic associated with the data analysis and 
the third is the systematic associated with the chiral extrapolation. 



6.2 AA-Scattering 

The simplest two-baryon system to write code for is AA, but this is also quite complex from a 
physics standpoint due to inelastic channels that become degenerate in the SU(3) limit. One should 
view this is an interpolating field for the / = 0, S = 2 channel. In fig. ^ we show the (preliminary) 
correlation functions we have obtained from some of the MILC lattices. The lightest quark mass 
correlator is quite noisy, and it will be difficult to extract a scattering amplitude from this data 
cleanly. It is clear that an increase in statistics by relatively small factors will allow for a extraction 
of the scattering lengths with uncertainties much less than ±1 fm at pion masses below ~ 500 MeV. 



7. Conclusions 

I am very excited by the fact that lattice QCD studies of nuclei and multi-hadron systems are 
an important part of the future of nuclear physics. However, lattice QCD calculations alone are 
not sufficient, and it is clear that a parallel development of lattice QCD and effective field theory 
efforts in this area is required for significant progress to be made. Despite the large length scales 
that characterize nuclear interactions, the presently available lattices allow for investigations into 
these systems, and we are attempting to do this presently. 

I would like to thank my collaborators for the many exciting discussions that generated this 
work. 
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Figure 9: Preliminary data: the logarithm of the ratio of the AA correlation function to the square of the A 
correlation function. (Each has been shifted vertically by an arbitrary amount for the sake of clarity.) 
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